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Abstract 

In this work we use the recently measured neutron occupancies in 
the 76 Ge and 76 Se nuclei as a guideline to define the neutron quasipar- 
ticle states in the lpOfOg shell. We define the proton quasiparticles by 
inspecting the odd-mass nuclei adjacent to 76 Ge and 76 Se. We insert 
the resulting quasiparticles in a proton-neutron quasiparticle random- 
phase approximation (pnQRPA) calculation of the nuclear matrix el- 
ement of the neutrinoless double beta (0vf3(3) decay of 76 Ge. A re- 
alistic model space and effective microscopic two-nucleon interactions 
are used. We include the nucleon-nucleon short-range correlations and 
other relevant corrections at the nucleon level. It is found that the re- 
sulting Oz//3/3 matrix element is smaller than in the previous pnQRPA 
calculations, and closer to the recently reported shell-model results. 
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The neutrinoless double beta (Ov/3/3) decay of atomic nuclei plays a key 
role in the search for massive Majorana neutrinos and their mass scale [1, 2, 3, 
4, 5, 6]. To extract the necessary information from the measured data [7, 8, 9] 
the nuclear-structure effects have to be accounted for by computation of the 
associated nuclear matrix elements (NME's). Here we discuss the NME's for 
the light-neutrino exchange mechanism. 

One most often uses the proton-neutron quasiparticle random-phase ap- 
proximation (pnQRPA) [3, 4, 10] to treat the structure of double-beta de- 
caying nuclei. Also recent shell-model results are available [11, 12, 13, 14]. 
The pnQRPA is a model tailored to describe the energy levels of odd-odd 
nuclei and their beta decays to the neighboring even-even nuclei [15]. Also 
its extension, renormalized pnQRPA [16], has been used [17, 18] to compute 
double-beta matrix elements. The pnQRPA (and the renormalized pnQRPA) 
contains a free parameter, the so-called particle-particle strength parameter, 
g pp , that controls the magnitude of the proton- neutron two-body interac- 
tions in the T = pairing channel [19, 20]. The fixing of the value of this 
parameter has been done either by using the data on two-neutrino double 
beta (2u(3/3) decay [18] or the data on single beta decay [21, 22]. 

The pnQRPA calculations are based on quasiparticle states [15] produced 
via a BCS calculation [23]. The BCS method gives occupations of the single- 
particle orbitals and the occupation amplitudes are connected to the energy 
differences between orbitals. The customary way to determine the single- 
particle energy difference is to use the Woods-Saxon mean-field potential 
[23]. Slight adjustments of the resulting energies can be done based on the 
data on energy levels of odd-mass nuclei in the neighbourhood of the nucleus 
where the pnQRPA calculation is done [24]. 

Recently the single-particle occupancies in the 76 Ge and 76 Se nuclei were 
measured by (p,t) reactions [25]. As a result, the vacancies in the neutron 
subspace lp-0f 5 / 2 -0g 9 /2, hereafter called the pfg subspace, could be deduced. 
In this work we exploit this spectroscopic data in pnQRPA calculations of 
the nuclear wave functions of the 76 As nucleus, the intermediate nucleus in 
the double beta decay of 76 Ge. We first produce the 2uPP NME to fix the 
allowed values of the g pp parameter by data on the half-life of the decay 1 . 
Based on this we finally compute the Ov(3(3 NME's. 

1 Since 76 Ge lacks single beta decay data we use the 2i/(3(3 data to fix the possible values 
of g PP . 
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Here we assume that the double beta decay of 76 Ge proceeds through the 
virtual states of the intermediate nucleus 76 As to the ground state of the final 
nucleus 76 Se. By assuming the neutrino-mass mechanism to be the dominant 
one, we can write the inverse of the half-life for the 0u(3(3 decay as [3] 

fgjr^Gf f^V(M^) 2 , (i) 

where m e is the electron mass and is the leptonic phase-space factor. 

The 0u(3(3 nuclear matrix element consists of the Gamow-Teller, Fermi 

and tensor parts as 

M<°"> = Mq$ - (^j 2 + • (2) 

Our numerical calculations verify the shell-model results of [14] and show 
that the tensor part in (2) is quite small and its contribution can be safely 
neglected in what follows. The expressions for the phase-space factor, the 
effective neutrino mass (m u ) and the matrix elements of (2) are given, e.g., 
in [1, 3, 10]. 

Table 1: WS single-particle energies for the pfg subspace in units of MeV. 





76 


Ge 


76 


Se 


Orbital 


neutrons 


protons 


neutrons 


protons 


iPs/2 


-11.52 


-9.015 


-12.77 


-7.000 


0f 5 /2 


-10.72 


-8.212 


-11.98 


-6.288 


!Pl/2 


-9.797 


-6.789 


-10.98 


-5.001 


0g 9 /2 


-7.030 


-5.311 


-8.305 


-3.371 



The nuclear-structure calculations are performed as described in [26, 27, 
28] . As a model space we have used the N = 3 and N = 4 oscillator shells and 
the 0hn/2 single-particle orbital, both for protons and neutrons. The starting 
values of the single-particle energies are obtained from the Coulomb-corrected 
Woods-Saxon potential, hereafter called WS, with the parametrization of 
[29]. These energies are presented for the pfg subspace in Table 1. The 
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measured neutron vacancies [25] in this sub-space have been summarized in 
Table 2. As can be seen from this table, the computed vacancies in the WS 
basis are far from the measured ones: the 0g 9 / 2 orbital is much too thinly 
occupied and the other orbitals are too full. 

The pnQRPA calculation of the wave functions of 76 As needs as input 
the neutron and proton occupation amplitudes in the pfg subspace. Instead 
of using the WS based BCS amplitudes for neutrons we now resort to the 
occupation amplitudes derived from the measured vacancies of Table 2. It is 
then interesting to see to what extent this choice of occupation amplitudes 
affects the magnitudes of the computed double beta decay NME's. 

We have to stress here that in this work we compute the occupancies 
of the neutron orbitals at the BCS level. Though effects beyond pairing 
may be present in the observed occupancies, in the present calculation we 
identify them with pair correlations, treated by the use of BCS. This can be 
justified partly by the good agreement between the calculated and observed 
properties of the low-lying energy levels in the adjacent odd-mass nuclei. 
Further support can be sought from calculations that determine the orbital 
occupancies and the pnQRPA amplitudes X and Y self-consistently [30]. In 
[30] this kind of scheme was applied in a realistic model space to study Fermi 
pn excitations in 76 Ge. From Fig. 5 of [30] one can see that the occupancies 
of the neutron orbitals in the pfg subspace do not change significantly when 
the BCS+pnQRPA results are replaced by the self-consistent results. The 
proton Fermi surface sharpens a bit when going from the BCS+pnQRPA over 
to the self-consistent approach of [30], but the effect is quite small (anyway, 
we do not have experimental occupancies available for protons). 

Table 2: Measured and WS based neutron vacancies in the pfg subspace. 





76 


Ge 


76 Se 


Orbital 


Exp. 


WS 


Exp. WS 


I/lpi/ 2 + ^lp 3 /2 


1.13 


0.357 


1.59 0.495 


^0f 5/2 


1.44 


0.500 


2.17 0.618 


^Ogg/a 


3.52 


5.43 


4.20 7.06 



No data on proton spectroscopic factors in the pfg subspace exist for Ge 
or 76 Se. The computed proton vacancies in the WS basis are given in Table 3. 
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One possible way to access the needed occupations for protons is to compute 
the spectra of the proton-odd nuclei 77 As and 77 Br, adjacent to 76 Ge and 
76 Se. For this one can use the quasiparticle-phonon coupling in diagonalizing 
the residual nuclear Hamiltonian. If one takes the phonons to be pnQRPA 
phonons one ends up with a scheme called the proton-neutron microscopic 
quasiparticle-phonon model (pnMQPM) [31]. The pnMQPM wave function 
consists of a linear combination of one-quasiparticle and three-quasiparticle 
components. A one-quasiparticle state in this scheme is a state whose wave 
function is dominated by a definite one-quasiparticle component. 

Table 3: Predicted proton vacancies in the pfg subspace. The calculations 
were done in the WS and adjusted basis. 





76 


Ge 


76 


Se 


Orbital 


Adj. 


WS 


Adj. 


WS 


!Pl/2 


1.82 


1.80 


1.70 


1.62 


!P 3 /2 


2.19 


1.92 


1.75 


1.44 


0f 5/2 


4.58 


4.43 


3.80 


3.32 


0g 9 /2 


9.28 


9.73 


8.61 


9.49 



Comparing the energies of the computed one-quasiparticle states with the 
available data gives information on the quality of the underlying BCS cal- 
culation and thus on the mean-field single-particle energies. The pnMQPM 
computed proton one-quasiparticle states in 77 As and 77 Br are given in Ta- 
ble 4. The wave functions corresponding to these states are indeed dominated 
by a single one-quasiparticle component. This is visible from the contribu- 
tion of this component to the total normalization of the wave function, as 
presented by the percentage in the last column of the table. 

The predicted one-quasiparticle energies of Table 4 correspond nicely to 
their measured counterparts. The critical proton orbit for the success of the 
pnMQPM calculation is vr0g 9 / 2 - The energy of this orbit had to be lowered 
to —7.1 MeV for 76 Ge and to —5.3 MeV for 76 Se to produce a reasonable one- 
quasiparticle spectrum. The corresponding energies we call adjusted energies 
and the computed proton vacancies in this basis are given in Table 3 in 
the columns called 'Adj.'. Comparison of tables 2 and 3 indicates that the 
adjusted proton and experimental neutron vacancies behave qualitatively 
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in the same way. However, the adjustments on the proton side have not 
such drastic effects on the vacancies as is the case on the neutron side when 
going from the WS to the experimental vacancies. For completeness, we also 
give in Table 4 the computed and measured energies of the neutron one- 
quasiparticle states in 77 Ge and 77 Se. Here the basis, adjusted independently 
of the pnMQPM calculation, shifts the 9/2+ state slightly below the measured 
1/2" state. Related to this one has to bear in mind that such very fine details 
of odd-mass spectra are hard to reproduce within a quasiparticle-phonon 
coupling scheme. 



Table 4: One-quasiparticle states in odd- A nuclei near 76 Ge and 76 Se. In the 
last column the contribution to the total normalization is given in per cents. 



Nucleus 


State 


E(exp.)[MeV] 


E(th.)[MeV] 


Main component (%) 


77 Ge 


1/2" 


0.160 


0.193 


z/lp 1/2 (94.7%) 




9/2+ 


0.225 


0.000 


z/0g 9/2 (97.0%) 




3/2" 


0-0.215 


0.014 


7rlp 3/2 (95.6%) 


77 As 


5/2" 


0.246 


0.000 


7r0f 5/2 (92.5%) 




9/2+ 


0.475 


0.471 


7r0g 9/2 (86.5%) 




1/2- 


0.000 


0.227 


z/lp 1/2 (93.6%) 


77 Se 


9/2+ 


0.175 


0.000 


z/0g 9/2 (97.6%) 




3/2" 


0.239 


0.677 


i/lp 3/2 (85.5%) 




5/2- 


0.250 


0.506 


i/Ofs/2 (93.0%) 




3/2" 


0.000 


0.083 


7rlp 3/2 (98.5%) 


77 Br 


9/2+ 


0.106 


0.072 


Tr0g 9/2 (93.3%) 




5/2" 


0.162 


0.000 


Tr0f 5 /2 (95.3%) 




1/2- 


0.167 


0.426 


7rlp 1/2 (87.6%) 



The main point of the present Letter is to see how the neutron vacancies 
of the pfg subspace, extracted from experimental data, affect the magnitude 
of the OvfiP NME's. As a side line one can also try to relate the measured 
vacancies to the single-particle energies of this subspace. The WS potential 
is a global parametrization of the nuclear mean field and is based on data 
on nuclei close to the beta stability line. This potential produces a smooth 
and gentle variation of the single-particle energies as a function of the proton 
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and neutron numbers. On the other hand, the WS potential is just an ap- 
proximate substitute for the Hartree-Fock (HF) or Hartree-Fock-Bogoliubov 
(HFB) methods to calculate the self-consistent mean field. The HF and HFB 
methods can produce results that are very different from those of the WS 
potential. Their results are also very much dependent on the two-body in- 
teraction used in the calculations (for recent articles on these features see 
[32, 33]). 

In [32] a discussion of the N = AO and Z = 40 mean-field gaps between 
the lp^ and 0g 9 / 2 single-particle orbitals was carried out by using the HF 
method with two different two-body interactions. Let us first discuss the 
situation for protons. It was found that on the proton side for one of the 
interactions, VMS (see Fig. 15 of [32]), the gap between the lp 1 / 2 an d 0g 9 / 2 
orbitals disappeares around Z = 34 and the two orbitals become inverted in 
energy. For the other interaction, JW (see Fig. 16 of [32]), the gap between 
these two orbitals diminishes when going from Z = 26 to Z = 48 but never 
disappears. Thus the VMS-computed mean field closely corresponds to our 
adjusted proton basis where the gap between the lp!/ 2 and 0g 9//2 orbitals 
disappears. 

On the neutron side the VMS based HF calculation from N = 24 to 
N = 48 produces a clear gap between the lpi/ 2 and 0g 9 / 2 orbitals, as seen 
in Fig. 3 of [32]. However, in [32] it was found that the effect of this gap 
disappeared when the T = 1, J = pairing interaction was added through the 
HFB method. If one would like to simulate these vacancies in the present 
simple WS+BCS calculations the WS mean field should be modified such 
that the gap at N = 40 closes. 

At this point it has to be stressed, however, that in the present work we 
do not have to resort to the details of the underlying mean field and the 
associated neutron single-particle energies since we have already the exper- 
imental vacancies available. Some further insight in the complexity of the 
self-consistent mean- field calculations for Ge isotopes is given in [33]. There 
several standard interactions were used in the Gogny-HFB and Skyrme- 
HF+BCS frameworks and very different results, from triaxial to axially sym- 
metric shapes, were obtained with different interactions for 76 Ge. In fact, 
looking at the simple Nilsson diagram indicates that a tiny oblate deforma- 
tion would suffice to close the N = 40 gap. 

After settling the problem with the occupation amplitudes of the single- 
particle states we are ready to compute the 2vf3f3 and 0vf3f3 NME's. As usual, 
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we consider the two extreme values of the axial- vector coupling constant, 
namely the bare value g\ = 1.25 and the strongly quenched value g& = 
1.00. When calculating the 0v{3(3 half-lives it is convenient to remove the 
dependence from the phase-space factor by redefining the NME as 

NME' = NME 

These redefined nuclear matrix elements are the ones that are listed in Tables 
5 and 6. 

In Table 5 we list the adopted g pp values as extracted by comparing the 
measured 2vf3f3 half-lives with the computed ones. By using these values of 
g pp we have calculated the Qv/3/3 NME's of (3) and we summarize their values 
in Table 5. In these calculations we have included the higher-order terms of 
nucleonic weak currents and the nucleon's finite-size corrections in the way 
described in [18, 34]. We have accounted for the short-range correlations by 
the Jastrow and UCOM (unitary correlation operator method) correlators, 
as discussed, e.g. in [26, 27, 28]. 



Table 5: Matrix elements of (3) computed in this work for different values of 
(?a an d g pp - For the short-range correlations both the Jastrow and UCOM 
prescriptions have been used. 









Jastrow 






UCOM 




9a 


9pp 








(Mg?)' 


(Mj? v) y 




1.25 


1.12 


2.288 


-0.772 


2.779 


3.385 


-1.143 


4.112 


1.00 


1.10 


1.700 


-0.579 


2.279 


2.413 


-0.818 


3.231 



Our computed results of Table 5 can be compared with the results of 
other recent works in the field. A selection of recent calculations including 
the Jastrow and the UCOM correlator is given in Table 6. The second and 
third columns of this table give the results of [28] where exactly the same 
methods as here were applied, the only difference being the use of a different 
set of single-particle energies, where the neutron Ogg/2 orbital was shifted a 
good one MeV to better reproduce the low-energy spectra of 77 Ge and 77 Se 
in a BCS calculation. By comparing the results of these two calculations in 
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Tables 5 and 6 one notices a significant reduction in the value of the total 
Ovfifi matrix element M^ u \ Furthermore, the Tubingen results [35, 36] are 
consistent with the results of [28]. 

Table 6: Values of the matrix element (M^)' of (3) obtained in some other 
recent works. Here (J) stands for Jastrow and (U) for UCOM. 



9A 


(J) [28] 


(U)[28] 


(J) [35, 36] 


(U)[36] 


(J) [14] 


(U)[14] 


1.25 


4.029 


5.355 


4.68 


5.73 


2.30 


2.81 


1.00 


3.249 


4.195 


3.33 


3.92 







The results of the shell model [14] are the smallest in Table 6. Interest- 
ingly, our present results for the Jastrow correlator, (M^)' = 2.779, and for 
the UCOM correlator, (M {0u) )' = 4.112, are closer to the shell-model result 
than the previous values quoted in [28]. The reduction of the magnitude of 
the pnQRPA calculated NME, which yields a value close to the shell-model 
result, is significant and deserves the further detailed study performed below. 

The reason for the reduction of the magnitude of the Qvf3f3 NME can 
be summarized by looking at the multipole decomposition of the NME. As 
an example we use the Jastrow correlated NME's. For the Fermi matrix 
element the reduction stems from the + intermediate states, as seen in 
Fig. 1. From Fig. 2 we see that for the Gamow-Teller matrix element Mq^ 
the significant changes concentrate on the 1 + and 2~ contributions. In the 
present calculations the 1 + and 2~ contributions are 0.448(1 + ) and 0.388(2") 
whereas in the previous work [28] they read 0.712(1 + ) and 0.899(2"). Thus 
the 1 + contribution has reduced by 37% and the 2~ contribution by 56%. In 
the old calculation the contribution of the 2~ was dominant but now it has 
reduced below the 1 + contribution. In both calculations the 1 + multipole 
has several important contributions whereas the 2~ contribution is coming 
almost solely from the first 2~ state. Hence the wave function of the 2f state 
plays a key role when seeking the reason for the reduction of the magnitude 
of the Gamow-Teller matrix element. 

The quality of the lowest 2~ state in the intermediate nucleus 76 As can be 
tested by computing the f3~ decay log ft values for transitions from this state 
to the ground state and one- and two-phonon states in 76 Se. The obtained 
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Figure 1: Multipole decomposition of the Fermi matrix element with Jastrow 
short-range correlations. The black bars correspond to the modified neutron 
and proton occupations. 



results are compared with the data and the calculations of [28] in Table 7. As 
can be seen there is a drastic improvement in the log ft value of the ground- 
state-to-ground-state transition. This transition tests exclusively the 2j~ wave 
function whereas the rest of the transitions depend also on the final-state 
wave function, built from the 2f collective phonon in the 76 Se nucleus. It is 
worth pointing out that in the present calculation the quasiparticle spectrum 
is more compressed than in the calculation of Ref. [28]. This increases the 
collectivity of the 2f state in 76 Se and thus results in smaller effective charges 
when trying to reproduce the data on the E2 transition probability from this 
state. 

The single f3~ decay is a non-trivial way to check the reduction in the 
0u(3(3 NME: The (3~ NME is reduced by 58% from the old value [28]. This is 
in nice agreement with the 56% reduction in the 2~ contribution to the Ovfifi 
NME. The main component of the wave function driving both transitions is 
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Figure 2: The same as Fig. 1 for the Gamow-Teller matrix element. 



the proton 0f 5 / 2 orbital coupled to the neutron 0g 9 / 2 orbital. In the present 
calculation the wave function of the 27/ state is more fragmented and thus 
reduces the pnQRPA amplitude responsible for the transitions. On the other 
hand, the occupation of the neutron 0gg/ 2 orbital has increased which also 
reduces the decay amplitudes since they are proportional to the emptiness of 
z/0g 9 / 2 . Similar considerations, though in a more complicated way, apply to 
the intermediate 1 + contribution. 

Our calculations show that the main contributions to the Qvf3f3 NME 
come from inside the pfg subspace. The implementation of the experimental 
occupations in the pnQRPA calculation brings the pnQRPA results closer to 
the shell-model results of [14]. The small contributions from outside the pfg 
subspace partly explain the deviations from the shell model result. In [37] 
the effect of expanding the shell-model single-particle basis was examined. 
Using the pfg subspace plus two-particle-two-hole excitations from the Ofy/ 2 
orbital it was concluded that these 2p-2h excitations increase the magnitude 
of by at most 20%. It still remains an open question how the differences 
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Table 7: Beta-minus decay log ft values for transitions from the 2 of 76 As 
to the ground state and one- and two-phonon states in 76 Se. 





o K + , 


n 


0+ 


9+ 


4f 


Exp. 


9.7 


8.1 


10.3 


8.2 


11.1 


Present calc. 


9.7 


7.4 


9.0 


8.4 


10.7 


[28] 


9.0 


7.7 


9.2 


8.7 


10.9 



between the shell-model and pnQRPA matrix elements tie to the omitted 
single-particle orbitals in the shell model and the shell-model occupations of 
the pfg subspace. 

Finally, our presently computed variations (M^)' = 2.279 — 2.779 (Jas- 
trow) and (M^)' = 3.231 - 4.112 (UCOM) in the 0u(3(3 NME can be 
converted to the following half-life limits 

4/? = (5.36 - 8.04) x 10 24 yr/((^)[eV]) 2 (Jastrow) , (4) 
4°/? = (2.45 -4.00) x 10 24 yr/((m^)[eV]) 2 (UCOM) . (5) 

In this Letter we have performed a pnQRPA calculation of the nuclear 
matrix elements involved in the neutrinoless double beta decay of 76 Ge. We 
have used a microscopic two-nucleon interaction in a realistic model space, 
and the calculations exploit the occupation amplitudes extracted from the 
recently available data on the neutron vacancies in 76 Ge and 76 Se. The sub- 
sequently calculated Qv/3/3 nuclear matrix elements are smaller in magnitude 
than the ones obtained in a standard calculation using the Woods-Saxon 
based single-particle occupations. This stems from the reduction in the con- 
tributions of the + , 1 + and 2~ intermediate states, with a special emphasis 
on the first 2~ state. These changes are related both to the revised occu- 
pations and to the changes in the pnQRPA amplitudes that derive from the 
revised occupations. 
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